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ABSTRACT 

We study the angular momentum evolution of binaries containing two white dwarfs 
which merge and become cool helium-rich supergiants. Our object is to compare pre- 
dicted rotation velocities with observations of highly evolved stars believed to have 
formed from such a merger, which include the RCrB and extreme helium stars. 

The principal case study involves a short-period binary containing a 0.6 Mq 
carbon-oxygen white dwarf, and a 0.3 Mq helium white dwarf. The initial condition 
for the angular momentum distribution is defined by the orbital configuration where 
the secondary fills its Roche Lobe. 

Since mass transfer from the secondary is unstable, the white dwarf breaks up 
on a dynamical timescale. After accreting some mass, the primary is assumed to 
ignite helium and evolve to become a yellow supergiant with a helium-rich surface. We 
assume conservation of angular momentum to compute the initial angular momentum 
distribution in a collisionless disk and subsequently in the giant envelope. At the end of 
shell-helium burning, the giant contracts to form a white dwarf. We derive the surface 
rotation velocity during this contraction. 

The calculation is repeated for a range of initial mass ratios, and also for the case 
of mergers between two helium white dwarfs; the latter will contract to the helium 
main-sequence rather than the white dwarf sequence. 

Assuming complete conservation of angular momentum, we predict acceptable an- 
gular rotation rates for cool giants and during the initial subsequent contraction. How- 
ever such stars will only survive spin-up to reach the white dwarf sequence (CO+Hc 
merger) if the initial mass ratio is close to unity. He+He merger products must lose 
angular momentum in order to reach the helium main sequence. 

Minimum observed rotation velocities in extreme helium stars are lower than our 
predictions by at least one half, indicating that CO+He mergers must lose at least one 
half of their angular momentum, possibly through a wind during shell-helium burning, 
but more likely from the disk following secondary disruption. 
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1 INTRODUCTION 

Following IWebbinkl Jl984) . ISaio fc Jeffervl (1200211 have 
demonstrated that the most probable origin for extreme he- 
lium stars is a stellar merger in a binary system containing 
a carbon-oxygen and a helium white dwarf, although some 
may originate in syste ms containing two helium white dwarfs 
JSaio fe Jeffervll200dh In the first case, the product ignites 
helium in a shell at the surface of the CO core and expands to 
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become a helium-rich supergiant. After the helium-burning 
shell burns outwards through most of the helium-rich sur- 
face layers, the star contracts to become a white dwarf. In 
the second case, the helium-ignition again occurs in a shell 
at the core-envelope boundary, but then burns inwards, lift- 
ing the electron-degeneracy in the helium core. When the 
helium-burning flame reaches the centre, the star essentially 
becomes a low-mass helium main-sequence star. 

A criticism occasionally heard (but, to our knowledge, 
not yet written) is that the products of such mergers should 
have a very high angular momentum and thus should be 
observed as rapid rotators. While general arguments suggest 
that this should not be a problem, it seemed appropriate to 
investigate the question more rigorously. 
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Previously ISaio fc Jeffervl (l200d . 120021) examined the 
evolution of the internal structure of stars following two 
types of white dwarf merger. In this paper we focus on the 
evolution of the angular momentum distribution from bi- 
nary progenitor to final white dwarf. We aim to set limits 
on the predicted rotation velocities of the merger products. 



2 ANGULAR MOMENTUM: GENERAL 
CONSIDERATIONS 

The final mass of systems that becom e extreme helium star s 
does not appear to exceed 0.9 Mq iSaio fc Jeffervl l2002ft . 
Assuming tidal locking for the progenitor binary we con- 
clude that the orbital angular momentum is dominant. We 
already know that angular momentum will be removed by 
gravitational radiation on a timescale that depends on the 
initial separation of th e binary and will lead to a decrease in 
the orbital separation llChaull 1978). For binary white dwarfs 
that have a period of a few hours and an initial separation 
of ~ 10 11 cm ~ 2R,0 the time needed to merge is less than 
the Hubble time. So we expect a significant fraction of exist- 
ing close white dwarf binaries to merge. Because of orbital 
decay, at some point the secondary fills its Roche Lobe. The 
orbital period at this point can be determined from the sep- 
aration and is approximately 3 minutes. Subsequently, mass 
starts to transfer from the secondary to the primary. The 
transferred mass cannot all be accreted onto the primary 
because of the Eddington limit. The mass transfer is un- 
stable because of the inverse mass-radius relation for white 
dwarfs. So the secondary will break up and form a thick 
disk round the primary on a dynamical timescale. For such 
systems, this is very short and of the order of the orbital 
period. 

To begin with, we assume conservation of total mass. 
This has been de monstrated for the early evolution of white 
dwarf mergers (jSe^retain^Cliabrier fc Mochkovitchl Il997t 
iGeurrero. Garcia^err^^^seTnll2^0^) /but not for the later 
stages of their evolution. In order to determine the char- 
acteristics of the angular momentum profile of the disk we 
use two quantitites, specific angular momentum and angular 
momentum density, defined in section 3. The conservation 
of total angular momentum is a plausible assumption. In 
addition, we assume conservation of angular momentum for 
every single particle (or fluid element). Thus the discussion 
reduces to an idealised collisionless model. 

This is a conservative assumption, which effectively re- 
duces to a worst case scenario. In the collisionless model, 
the transfer of angular momentum is slow compared to the 
dynamical timescale, so, it does not have important effects. 
The opposite assumes a model in which collisions are im- 
portant, producing a viscous disk. The problem of angular 
moment um transport in a viscous disk has been approached 
already (iLvnden-Bell fc Pringlelll974l) . where it was shown 
that "the angular momentum is steadily concentrated onto 
a small fraction of the mass which orbits at greater and 
greater radii while the rest is accreted onto the central body" 
on an approximately dynamical timescale. The problem for 
the viscous disk approximation is that the accretion rate for 
the surviving white dwarf then becomes very high and it is 
not clear whether subsequent evolution will lead to the as- 
similation of this mass or its complete ejection by radiation 



pressure. Consequently, we have started with the collision- 
less model to see whether it, too, faces major obstacles. 

Using the above assumption we propose an angular mo- 
mentum distribution profile for the disk. In the purely colli- 
sionless model, as will be shown, the inner radius of the disk 
is approximately twice the radius of the primary. The inter- 
vening gap will prevent accretion. In real systems, random 
collisions will broaden the angular momentum profile. We 
will simulate this by applying a Gaussian smoothing func- 
tion, where the width of the Gaussian may be considered 
representative of the frequency of the collisions. Note that 
this procedure does not change the total angular momentum 
of the system. 

Since the disk was created from the decomposition of 
the secondary, which was a helium white dwarf, helium 
will be its dominant constituent. Following Saio fc Jeffervl 
(2002), after a small quantity (0.004 Mq) of helium is ac- 
creted by the primary, helium ignition occurs at the base 
of the accreted layer. This energy source forces the star to 
expand to become a giant. Models indicate that the radius 
of the giant will be two orders of magnitude greater than 
that of the disk, initially about 60 Rq. 

A simple model for the angular momentum distribution 
in the giant can be obtained by assuming that each cylindri- 
cal element in the disk forms a spherical shell conserving its 
angular momentum. In order to describe the density profile 
of the giant we assume it consists of a degenerate core, that 
is the carbon-oxygen white dwarf, and a convective envelope. 
In the fully conservative case, the final mass of the convec- 
tive envelope is equal to the mass of the initial disk, and 
hence of the helium secondary. Therefore we can describe 
an angular velocity distribution for the star. 

Following calculation of the angular velocity distribu- 
tion in the giant merger product, we investigate how contrac- 
tion affects the rotation. After helium-burning is completed, 
the stars will contract, and hence rotate more quickly. The 
question is whether their rotation will approach the critical 
breakup velocity and what the consequences might be. 

For the CO+He WD merger, we investigate three pos- 
sible cases. In the first case, the central region rotates as a 
rigid body and the angular velocity profile of the envelope 
depends on the initial conditions, meaning the disk angular 
momentum distribution. In the second case, 25% of the con- 
vection zone near the surface rotates as a rigid body as well 
as the central region whereas there is differential rotation 
in the intermediate region. Thirdly, we examine the case of 
completely rigid body rotation. This should be the ultimate 
equilibrium state since no shear to rques occur that could 
lead t o angular mo mentum transfer iLvnden-Bell fc Pringld 
Il974t |Pringi3ll98lT) . however it is not likely to be achieved, 
since it is a slow process. 

Using the above assumptions we will model the an- 
gular momentum evolution of a binary which consists 
initially of two white dwarfs. For the structure of the 
white dwarfs, we w ill adopt the models described by 
IChand rasekhari il95STl . For the principal calculation we will 
consider the primary to be a carbon oxygen (CO) white 
dwarf with mass 0.6 Mq and radius 0.013-Rq and the 
secondary a helium ( He) white dwarf with mass 0.3 
and radius 0.021 |V ennes 1 F ontaine, fc Brassardlll995l : 
iPanei. Althaus fc Benvenutol200ol) . In addition we will con- 
sider CO+He binaries with masses of 0.7Mq + 0.2Mq and 
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0.5Mq+0AMq, in order to see the dependence, if any, of the 
angular velocity on the initial mass ratio. We also compute 
appropriate quantities for a number of He+He white dwarf 
configurations. Finally, the results of these calculations are 
compared with observed angular velocities in extreme he- 
lium stars. 



3 ORBITAL DECAY TO ROCHE LOBE 

According to the General Theory of Relativity, two orbiting 
masses Mi, M2 with a separation a will rad iate angular 
momentum at a rate jLandau fc Lifshitzlll958() 



where q is the ratio of the mass of the primary over the 
mass of the secondary. 

For the 0.6 + 0.3Mq system defined above, the sec- 
ondary will fill its Roche lobe when the orbital separation 
is 0.067 Rq. From eqns. 6 and 1 we can estimate the or- 
bital decay timescale. Assuming an initial separation of 1.5 
Rq, which corresponds to an orbital period of 5 hours, we 
obtain a timescale of 4.10 9 y, which is less than a Hubble 
time. We note that observations o f an increasing number 
of such close binary white d warfs jPaczvnskil Il990| : iMarshl 
Il995l iNapiwotzki et al.ll2005h are commensurate with esti- 
mated l^neTg^r^ates~Jlbe^^^i]^i3i iNelemans et al]l200ll) 
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where G is the gravitational constant, c the speed of 
light, and M = Mi + M2 is the total mass of the system. 
Due to angular momentum loss their orbits will decay. The 
total angular momentum of a tidally locked system will be 



Jtot — J{ + Ji + J2 + Ji 



(2) 



where the superscripts o and s refer to orbital and spin 
angular momentum respectively and the subscripts 1 and 2 
refer to the primary and to the secondary. We can express 
the orbital angular momentum in the frame of reference of 
the centre of mass 
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where M is the total mass of the system. Spin angular 
momentum can be expressed 
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where u> is the angular velocity and I refers to the mo- 
mentum of inertia. Moments of inertia can be evaluated from 
the density profile of the stars . However some distortion s 
may occur that alter it slightly llameslll964l : lTassou]|ll97Sl) . 
The angular velocity is given by 



GM 



(5) 



Substituting angular momentum expressions from equa- 
tions 3, 4 and differentiating with respect to a, we obtain 
the following expression 
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As the orbit decays, the equipotential surface surround- 
ing the two stars (Roche lobes) shrinks until one component 
exactly fills its own lobe. The Roche lobe radius is given by 
feggletonlll9"83l) 



OiL = 
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4 DISK FORMATION 

When the secondary fills its Roche Lobe it disintegrates and 
its remnants form a disk. We will use the assumption of an- 
gular momentum conservation to determine the mass distri- 
bution in the disk. First we construct two useful quantities 
for the distribution of angular momentum in the secondary. 
The first is the angular momentum per unit mass, here- 
after specific angular momentum, which will be expressed 
as dJ(m)/dm, assuming the system shows cylindrical sym- 
metry. Since we generally assume a collisionless model, this 
quantity remains constant throughout the problem unless 
some mass is lost. The other quantity is the angular mo- 
mentum contained within a distance r from the centre of 
the system, hereafter angular momentum density, and will 
be expressed as dJ(r)/dr. It depends strongly on the ge- 
ometry of the problem. When the secondary fills its Roche 
Lobes, it expresses the angular momentum contained in a 
thin slice of the star at distance r from the rotation axis of 
the primary. The slice lies perpendicular to the radial vector 
that points from the centre of the primary to the centre of 
the secondary. In the case of the disk, dj/dr expresses the 
angular momentum carried by a thin cylindrical shell of ra- 
dius r. In the case of the giant star, it expresses the angular 
momentum carried by a spherical shell at distance r from 
the rotation axis of the star. 

From the definition of angular momentum a particle of 
mass dm that lies at distance r from the primary and moves 
at angular velocity ui has angular momentum 



dj = r udra 



(8) 



Integrating the specific angular momentum over mass 
and the angular momentum density over distance we ob- 
tain the total angular momentum of the system. For our 
example conditions it is 4.5 x 10 50 ergs and remains con- 
stant throughout the evolution of the system, unless some 
mass is lost. We can easily obtain from the orbital elements 
the angular velocity of the secondary when a = oll to be 
uj — 0.0346 rads -1 . The density of the secondary can be 
approxi mated by a polynomial function that fits the data 
given by IChandrasekhail (ll958T> for the structure of white 
dwarfs. Thence we can derive the specific angular momen- 
tum (Fig.0 and the angular momentum density (Fig.|5J for 
the system immediately before disruption of the secondary. 

After disruption of the secondary, a disk is formed. Since 
the disk is supported by gravity we assume it to be Keple- 
rian. Therefore the angular velocity distribution can easily 
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Figure 1. The specific angular momentum dj/dm. It remains 
constant throughout the problem, since we examine a collisionless 
model. 
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Figure 2. The angular momentum density dj/dr at the Roche 
Lobe stage. We can see that most angular momentum is carried 
at a distance slightly further from the centre of the secondary 
(O.O67i?0). This is because dj/dr is proportional to the mass 
contained inside a slice of the star and to the square of the dis- 
tance. This product has a maximum at 0.069-Rq from the centre 
of the primary. 



be determined to first order, ignoring the self gravitating 
effects of the disk itself. 

/ ^ l GMl f«\ 
= ( 9 ) 

The angular momentum for a particle of mass m per- 
forming circular motion of radius r, is J = mur . Parti- 
cles that carry less of the angular momentum at the binary 
stage lie at the front surface of the star {i.e. face-on to the 
the primary). The angular momentum of these particles will 
determine the inner radius of the disk. In our case, this is 
rin = O.O22/?0. We can use the same argument for the outer 
radius of the disk. In that case we evaluate the angular mo- 
mentum of a particle at the rear surface of the secondary. 
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Figure 3. The angular momentum density dj/dr of the disk. An- 
gular momentum is spread in a more extended region than it was 
at the Roche Lobe stage. If we take into account the self grav- 
itation of the disk then its size decreases, but the total angular 
momentum cannot change. Since we have a collisionless model, 
no mass is scattered closer to the star than n n = 0.22Rq. The 
spin angular momentum of the primary has been neglected since 
it is three orders of magnitude smaller than the angular momen- 
tum of the disk. The dotted line corresponds to the convolution 
of dj/dr with a Gaussian with a = 0.005Rq/\/2. 



We conclude that it is r out = 0.305i?Q. If we take into ac- 
count the self gravitating effect of the mass stored in the 
disk, the outer radius is decreased 30% to r out = 0.203i?Q, 
whereas the inner radius will be unaffected. This may be 
shown by a first-order approximation. Assume that a parti- 
cle of unit mass carrying angular momentum J lies in the 
gravitational field of a mass M. At radius r the attractive 
force is GM/r 2 . Including the gravitational field of a disk of 
mass m, the central force on the test particle, now at radius 
r', is G(M + m)/r' 2 . Since J is constant, the ratio of radii 
(for m = O.3M , M = O.6M ) is r'/r = 2/3. In higher 
order, the geometry of the disk must be considered, but this 
does not affect the angular momentum calculation. 

Whether we use the first or the second approximation 
for the angular momentum density will have no effect on 
our final results for the rotation of the giant star, since 
the specific angular momentum does not change. The choice 
only affects the angular momentum density at the disk stage 
(Fig. E|. 

During the formation of the disk collisions may take 
place. These will broaden the angular momentum distribu- 
tion, in particular to scatter material into the region r < n n - 
Angular momentum will also be transfered to outer parts 
of the disk which will become more extended. In order to 
simulate the effect of these collisions we convolve the angu- 
lar momentum functions (dj/dr) with a Gaussian of width 
a = 0.005Rq /V3) where the width of the Gaussian may 
be taken to represent the overall efficiency of the collisions. 
The resulting density distribution demonstrates that only a 
very modest collisional redistribution of angular momentum 
is required to bring disk and star into contact (Fig. [3J . 
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5 GIANT STAGE 

Following disk formation and providing gas has been scat- 
tered to r < ri n , helium will be accreted from the disk onto 
the surface of the former pri mary, which becomes t he degen- 
erate core of the merged star. lSaio fc Je ffcrv (2002) adopted 
an accretion rate of roughly half the Eddington rate. After 
~ 0.029Mq is accreted, a helium flash occurs. This leads 
the star to expand and form a yellow giant. The expansion 
is rather rapid and lasts only 200 years. The radius of the 
giant initially reaches ~ 60-Rq . This stage of evolution may 
correspond to some hydroge n-deficient carbon gia nts and the 
coolest of the RCrB stars JSaio fc Jeffervll2002ft . The radii 
of some models do exceed this value, so we surmise that the 
surface rotation will be lower than that deduced here. How- 
ever, assuming a conservative distribution of specific angular 
momentum, spindown during expansion and spinup during 
contraction should lead to a similar final result. 

Our objective is to predict the angular velocity profile 
for such a giant. In order to model its density profile, we 
expect it to consist of a deg enerate core and a con vective 
envelope. Using tables from Ichandrasekharl il95Sf) . fitted 
with a 5 th order poynomial, we approximate the envelope 
as a polytrope of index n = 3/2 and express the density as 
a function of r that fits the nu merical results of the Lane - 
Emden equation for n = 3/2 iMohan fc Al-Bavatvlll98Cft . 
We assume that initially no angular momentum is trans- 
ferred between the various parts of the star. Therefore we 
first make the approximation that the outside surface of the 
star will form the outside shell of the star, whereas the inner 
parts of the disk will stay near the core, and that there is 
no change in the specific angular momentum. 

However, it is also possible that helium-shell ignition 
occurs before the disk has been fully accreted, so that the 
material which was accreted first becomes the surface of the 
expanding giant, while the disk survives for sometime inside 
the giant envelope. Subsequent accretion, or the disintegra- 
tion of the disk, feeds mass into the interior of the giant, 
close to the core-envelope boundary. We will examine this 
case as well, assuming that the outer layer of the giant, with 
a mass equal to the burning shell, comes from the inner part 
of the disk and the remainder of the envelope comes from the 
disentegration of the disk and obeys cylindrical symmetry. 

We neglect differential rotation within a shell, so a thin 
spherical shell of mass dm at distance r from the centre of 
the star rotates as a rigid body with momentum of inertia 



I = ^r 2 dm (10) 



The choice of cylindrical shells, rather than spherical, 
is made primarily on grounds of mathematical simplicity. 
There is some evidence (e.g. the Sun) that rotation may be 
a function of latitude - possibly suggesting cylindrical sym- 
metry. It will be seen, however, that the difference between 
the extreme cases of rigid-body and differential rotation are 
sufficiently small that the choice of geometry in the latter 
are unlikely to be important. 

We examine a number of cases. 
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Figure 4. The angular momentum density dj/dr of the convec- 
tive envelope of the giant. 
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Figure 5. The ratio of centrifugal force to gravitation force as a 
function of radius for case 1. The plot shows only the differentially 
rotating stellar envelope. 



Case 1: Rigid core + differential envelope 

The first case involves two regions including a central re- 
gion containing the degenerate core and a shell of the star's 
envelope of width 0.045 Rq performing rigid body rotation. 
Assuming that the core rotates with the original angular ve- 
locity of the primary white dwarf, this is the critical radius 
for rigid-body rotation. Beyond this, the envelope is assumed 
to rotate differentially according to the angular momentum 
distribution determined by the disk stage (Fig.[IJ. The outer 
boundary to the region of solid body rotation is chosen in 
order to avoid the possibility of having a rotational velocity 
that exceeds the orbital velocity at this point. This scenario 
predicts an angular velocity on the outer layer of the star of 
1.15 x 10 _7 rads _1 , corresponding to an equatorial surface 
velocity v cq = 4.95 km s -1 . There is a slight discontinuity 
between the two regions described above. It can be consid- 
ered as a surface of infinite gradient du/dr. In such a case 
the transfer of angular momentum will be very fast and the 
discontinuity will vanish. This will be the initial stage in the 
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way between the extremes discussed in cases 1 and 3. We ex- 
pect that near the surface, convection may be strong enough 
that mixing will cause these outer layers to rotate as a rigid 
body (Fig.JTJ. Therefore the angular velocity for the surface 
and the outer convection zone is 1.27 x 10 _7 rads _1 , corre- 
sponding to v cq = 5.33 kms -1 . There is also a slight discon- 
tinuity at the 45 Rq boundary of 5ui — 1.2 x 10 -8 rads -1 . 
This represents an intermediate case, since angular momen- 
tum is transferred within the star. 
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Figure 6. Logarithmic plot of the angular velocity ui at the giant 
stage as a function of the distance from the rotation axis. The 
central region rotates as a rigid body, whereas angular velocity 
falls with distance from the axis. 
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Case 3: Rigid body rotation 

The third case to consider is that of rigid body rotation for 
the whole star. Shear torques may transfer angular momen- 
tum outwards from the core to the outer layers and thus 
make them move faster, while the inner layers will rotate 
more slowly. Rigid body rotation is the equilibrium state 
since it is the lowest energy state for an object with de- 
termined ang ular momentum llLvnden-Bell fc Pringldfl974l : 
|Pringlelll98ll i and contains no shear torques that will trans- 
fer angular momentum. However, this equilibrium is not 
achieved by real stars even on very long timescales, as il- 
lustrated by the Sun. It is unlikely to occur on the very 
short timescales and very low-density envelopes under con- 
sideration here. Nevertheless, it represents an important lim- 
iting case in which most of the angular momentum is car- 
ried in the outer layers of the star. We can easily evalu- 
ate the moment of inertia for the giant star by integrat- 
ing Eqn. 10 throughout the star to obtain /giant = 2.85 x 
10 2 gcm 2 Mq -1 Rq -2 . From conservation of angular mo- 
mentum we find that J g i ant = 4.71 x 10 -5 ergs M Q -1 Rq -2 . 
Therefore 
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Figure 7. The angular velocity ui of the disk, in the second sce- 
nario. The central and the intermediate region rotate as in the 
previous case (Fig. |HJ, however, there is a layer at the outter part 
of the convective envelope that rotates as a rigid body. There are 
two slight dicontinuities at the boundaries of the three regions. 



giant's angular momentum evolution (FigEJ- Fig. shows 
the ratio of centrifugal force to gravitation force as a function 
of radius. This drops rapidly outside the core, so rotatation 
should not significantly affect the structure of the envelope, 
although it may produce a slightly ellipsoidal star. 

Case 2: Rigid core + transition layer + rigid 
envelope 

In the second case we divide the star into three regions, 
the first is the core and a thin layer that rotates as a rigid 
body as before. The second region lies between 0.045 Rq and 
45 Rq, where the star rotates differentially as before. The 
third region lies between 45 Rq and the surface. This bound- 
ary was chosen arbitrarily in order to represent a case some 



Jgiant 

which is 1.64 x 10 -7 rads -1 corresponding to v cq — 
7.1 kms -1 . This case represents the fastest rotation of the 
surface layers and sets an upper limit for the observed rota- 
tion velocity. 

We conclude that the equatorial surface velocities do 
not depend strongly on the details of the model. Using the 
above argument we can define an upper limit on the an- 
gular momentum that corresponds to rigid body rotation. 
Any velocity observed should be less than this limit. We 
have not taken into account any losses of angular momen- 
tum due to mass ejection. Such phenomena will lead to lower 
angular momentum and therefore lower angular velocity. In 
the polytropic density profile we have neglected any term 
from rotation; this gives a good approximation since the 
linear velocity of circular motion due to the gravitational 
field of the star at a radius equal to the giant star radius 
is 53 kms -1 , or about one order of magnitude greater than 
the velocities we found. In addition, any distortion caused 
by rotation will increase the moment of inertia of the star 
and decrease the angular velocity, supporting the statement 
that rigid body rotation of a spherical body represents the 
maximum possible equatorial surface velocity. 
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Figure 8. Logarithmic plot of the angular velocity lj as the star 
contracts by a fraction j3 = Ro/R- Ro is the radius of the giant 
before contraction. 



Case 4: Rigid core O.87M0 + differential envelope 

The fourth case can be considered as a combination of cases 
1 and 3. During evolution as a giant, the helium-burning 
shell eats up the envelope, so the core mass increases, while 
the core radius remains virtually unchanged. Conversely, the 
envelope retains the same radius but a drastically reduced 
mass. We consider a final configuration in which the enve- 
lope mass is O.O3M0. The moment of inertia of the core 
increases but cannot exceed 1.4 x 10 -5 gcm 2 . Even in this 
case, no more than 15% of the total angular momentum can 
be stored in the core before the core reaches its breakup 
limit. Therefore the envelope will have to rotate about 30 
times faster than in case 1. 

Such rotation is extremely fast and seems implausible; 
at least, it is not observed. The problem is that, to be conser- 
vative, angular momentum must be transferred into the core 
at the same time as the envelope is ingested. However, once 
the core reaches breakup velocity the transfer of angular mo- 
mentum becomes impossible, so the momentum must stay 
in the envelope, creating a paradox. Interestingly, this par- 
allels the case for the viscous disk IILvnden-Bell fc Pringla 
I1974T) in which angular momentum is expelled outwards as 
mass migrates inwards. 

This calculation therefore strongly suggests that if sta- 
ble shell-burning giants are to be produced by white dwarf 
mergers, then a large fraction of the total angular momen- 
tum must be dissipated from the disk before it is accreted 
onto the primary. 

While this case is the most interesting and realistic of 
the four considered, it has presented a paradox for the trans- 
fer of angular momentum from the envelope to the growing 
core. The consequences of this deserve more attention than 
we are currently equipped to explore. 



6 HOMOLOGOUS CONTRACTION 

Having considered the angular momentum distribution of 
a star as it evolves from a white dwarf binary through to 
being a giant, we next consider the angular velocity distri- 



Figure 9. The ratio of the velocity due to rotation over the 
velocity of the circular orbit on the surface of the star A, as a 
function of /3 which is the contraction factor for three pairs of 
masses of the progenitors, from top to bottom 0.7Mq + 0.2Mq, 
0.6Mq +O.3M0 and 0.5M & +0AM & . Contraction cannot con- 
tinue if A approaches unity. In order to become a white dwarf the 
star needs to contract by a factor of about 5.10 3 . We can see that 
for progenitor masses of 0.5Mq + 0.4Mq, contraction can reach 
the white dwarf radius, whereas for the other two cases angular 
momentum needs to be removed. The dotted line corresponds to 
the model of a merger with progenitor masses of 0.6Mq + 0.3Mq 
where 50% of the secondary mass is lost. 



bution following shell helium extinction and the subsequent 
contraction at constant luminosity towards the white dwarf 
cooling track. 

Since the core is degenerate, we need only consider the 
evolution of the envelope, which we shall consider to con- 
tract homologously. A simple model can be developed by 
assuming that the star rotates as a rigid body (case 3). 

However, no crucial differences will be found if we as- 
sume differential rotation, as demonstrated in the previous 
section. The only effect would be to find lower surface equa- 
torial velocities. In order to become a white dwarf the star 
needs to reduce in size by a factor ~2x 10 -4 . Leaving the 
core unaffected, while the star shrinks, we evaluate the angu- 
lar velocity and the equatorial surface velocity corresponding 
to this rotation. We plot the angular velocity (Fig. |HJ and 
the ratio of the linear equatorial rotation velocity w cq to the 
linear velocity of a particle moving on a circular orbit u or b 
at the surface: 

A = Ucq/Worb (12) 

(Fig. El middle curve) . The quantity on the horizontal axis is 
the ratio j3 of the radius of the giant Ro to the radius during 
contraction R. The momentum of inertia decreases with the 
square of the radius, which will lead the star to rotate faster 
(Fig. 0. Note that contraction implies increasing (3. 

Problems occur when the rotation velocity approaches 
or exceeds the value of velocity for a circular orbit (A > 1). 
Thus, when the star contracts to ~ 0.007 times its initial 
radius, the two velocities become equal and the star cannot 
shrink anymore unless some angular momentum is lost. At 
this point the velocity is of order ~ 10 2 kms -1 . This is much 
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Table 1. Theoretical predictions for the rotation of helium main-sequence stars produced by He+He white dwarf mergers for various 
combinations of initial masses. 



M/Mq 


Mi + M 2 /M Q 


//io- 3 m r2 ) 


J/lO-^rgsM^R" 2 


oj/rad s 1 


t>rot/10 3 kms 1 


f O rb/10 3 kms 1 


0.5 


0.3+0.2 


1.16 


2.53 


0.044 


3.2 


0.95 


0.5 


0.4+0.1 


1.16 


1.71 


0.030 


2.2 


0.95 


0.7 


0.4+0.3 


3.18 


3.83 


0.012 


1.2 


0.96 


0.7 


0.5+0.2 


3.18 


3.18 


0.010 


1.0 


0.96 



lower than the speed of light and relativistic effects may be 
ignored. 



7 AN ANGULAR MOMENTUM PROBLEM? 

So far we have examined a fully conservative scenario from 
the Roche lobe stage and after for a 0.6 + 0.3Mq binary. 
However, much energy is ejected at the breakup of the sec- 
ondary. Some of this energy may be deposited in the disk 
as thermal kinetic energy, while some material of high an- 
gular momentum is likely to leave the system, although the 
maximum possible ma ss lost has been shown to be small 
<Han fc Webbinldll99Sft . 

We therefore examine whether the contraction of the 
giant can proceed conservatively if some mass (and hence its 
associated specific angular momentum) is ejected. Although 
this is contrary t o our initial assumptions and previous nu- 
merical results jSesxetsM^Gh^rieL^Mod^ovit^ [iQQt]; 
lHan fe Webbinkl Tl^99T ~ lG^urrero^Glaxia^er^^gbernl 
2004), this mass ejection is treated as taking place during 
the phase of disk formation. This approximation avoids 
numerical complications which are beyond the scope of 
this paper, but allows us to obtain a rough estimate of 
the resulting configurations when mass loss is taken into 
account. 

We have solved the problem for the cases where 1%, 5%, 
10%, 25% and 50% of the mass of the secondary is ejected 
during the merger process. In all cases, the giant cannot 
contract to the radius of a white dwarf. Fig. includes the 
case in which 50% of the mass of the secondary is ejected. 
In this test case, 2/3 of the total angular momentum is re- 
moved, since particles carrying more angular momentum are 
more likely to escape. A reaches unity when (3 « 1000, which 
means that the star can reach a small size (10~ 3 times its 
original), but not less than five times the required white 
dwarf radius. 

In order to investigate the role of the initial conditions, 
we have solved the problem for various relative masses of 
the white dwarf binary components. We have plotted the 
ratio of the equatorial rotation velocity to the surface or- 
bital velocity (A) for progenitor binaries of 0.7 + 0.2Mq and 
0.5 + 0.4M Q in addition to the test case of 0.6 + 0.3Mq 
(Fig. EJl. When A > 1, contraction cannot proceed. We see 
that for mass ratios nearer unity the star can contract to 
white dwarf size without the need for additional angular 
momentum loss. This is expected, since such systems need 
to lose more angular momentum before the secondary star 
fills its Roche Lobe. This is because the white dwarf radius 



varies inversely with mass, so more massive white dwarfs 
will have smaller radii and will fill their Roche lobes later. 
Close orbits involve less angular momentum of course, and 
their products will rotate more slowly. 

Therefore conservation of angular momentum does 
present a problem for the white dwarf merger model in most, 
but not, all cases. The problem arises either during helium- 
shell burning, when angular momentum has to be trans- 
ferred to the envelope to avoid the core reaching breakup 
velocity, or during contraction, when the surface reaches 
breakup velocity. 



8 HE+HE WHITE DWARF MERGERS 

The case of merger between two helium whit e dwarfs is ini- 
tially identical to that of the CO+He merger iSedo fe Jeffervl 
2000). After shell helium ignition, the star expands to be- 
come a giant, and the envelope is capable of storing sub- 
stantial angular momentum. However, the star evolves to 
become a h elium main-sequence star on a short timescale 
(Ibcn 1990). The case is analogous to the CO+He merger 
which must spin up during contraction to become a white 
dwarf. 

We have calculated rotation rates for helium main- 
sequence stars arising from a number of initial configura- 
tions. The final configurati on is an n = 5/2 polytrope with 
overall dimensions given bv lPaczvnskil lll97ll) . 

TableQshows the run of the moments of inertia 7, total 
angular momenta J, and angular and linear equatorial ro- 
tation velocities (u},v to t) for a series of inital configurations 
(Mi, M2). The predicted values for v ro t are very high and, in 
most cases, substantially exceed the critical orbital velocity 
at the surface (v rb). Clearly such stars cannot exist and it 
is necessary to find a mechanism for expelling angular mo- 
mentum at an earlier stage. In this case losing some angular 
momentum at the breakup of the secondary when it fills its 
Roche Lobe is likely, since the total mass of the binary and 
the binding energy is lower. In addition to this, less massive 
white dwarfs have larger radii and fill their Roche lobe ear- 
lier. This will enable mass carrying angular momentum to 
escape the system earlier. 



9 OBSERVED ROTATION RATES 

Observations suggest modest rotation velocities for extreme 
helium stars. Measurements of the projected rotation veloc- 
ity w eq sin i for ten objects have been obtained by fitting line 
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Table 2. Projected rotation velocities for extreme helium stars. 
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Star T c ff/K logg i) eq sini/kms 1 /J/Rq oj sin i/10 6 rad s 1 Reference 



FQ Aqr 


8750 ± 300 


0.3 ±0.3 


20 ± 5 


136 ± 96 


0.22 ±0.17 


Pandev et al. (20061 


HD168476 


13500 ± 500 


1.6 ±0.25 


25 ± 5 


23 ± 13 


1.6 ± 1.0 


Pandev et al. (20061 


LSS 99 


15330 ± 500 


1.9 ±0.25 


30 ± 5 


16 ±9 


2.9 ± 1.7 


Jefferv et al. (19981 


HD124448 


15500 ± 500 


1.9 ±0.25 


4±5 


16 ±9 


0.38 ±0.5 


Pandev et al. (20061 


LSS 4357 


16130 ± 500 


2.0 ±0.25 


45 ± 5 


14 ±8 


4.9 ± 2.8 


Jefferv et al. (19981 


LS 11+33° 5 


16180 ± 500 


2.0 ±0.25 


45 ± 5 


14 ±8 


4.9 ± 2.8 


Jefferv et al. (19981 


V1920 Cyg 


16300 ± 900 


1.7 ± 0.25 


40 ± 5 


23 ± 13 


2.6 ± 1.5 


Pandev et al. (20061 


BD+10°2179 


16900 ± 500 


2.55 ± 0.2 


18 ±5 


6±3 


4.7 ± 2.5 


Pandev et al. (20061 


LSE 78 


18000 ± 700 


2.0 ±0.1 


20 ± 5 


15 ±3 


2.0 ±0.7 


Jefferv (19931 


DY Cen 


19500 ± 500 


2.15 ± 0.1 


20 ± 5 


13 ±3 


2.4 ±0.8 


Heber (19931 
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Figure 10. Projected angular rotation logu sin i as a function of 
relative inverse radius /3 for extreme helium stars. The dotted line 
represents a homologous contraction, i.e. u oc 1/R 2 , and the solid 
line represents the maximum theoretical rotation rate shown in 
Fig.! 



profiles to absorption lines in the optical spectrum (Tabled . 
For small values of u cq sini, these may been overestimated 
because of the difficulty of deconvolving instrumental and 
rotation broadening profiles in the observed spectra. Ta- 
ble H also gives the measured values of T e g and logg. By 
assuming a core-mas s shell-luminosity (M c — L s ) relation 
JSaio fc Jefferv! Il988l cf.), the latter may be combined to 
estimate both a mass and a radius for the extreme helium 
stars. Hence we can obtain the projected (minimum) angular 
rotation velocity ojsini. We can normalise the relative radii 
by assuming a minimum T e ff ,o before contraction at constant 
luminosity and hence derive /3 = Ro/R = (r o ff/T , e ft,o) 2 - 
We have adopted T c ff ,0 = 5 000K, the value at which post- 
merg er giants commence their contraction iSaio fe Jefferv! 
2002). A different value would impose a horizontal offset 
onto Fig. 1101 which demonstrates an overall increase in 
ajsini with increasing /3 (and hence with increasing T e ff)- It 
is interesting that the maximum slope represented by these 
data corresponds to the spin-up of a homologously contract- 
ing star, lu oc R~ 2 , as theory predicts. 




0.2 0.4 0.6 0.! 

log/3 



1.2 1.4 



The projected rotation rates are all smaller than the 
maximum rate predicted by theory assuming complete con- 
servation of angular momentum through the white dwarf 
merger. Assuming a random distribution of inclination an- 
gles, the average rotation rates should be n/4 times the max- 
imum rates. The observed rates are less than one third of the 
maximum rate. This implies that at least half of the angular 
momentum must have been lost, either during the merger 
process or whilst the star was a cool giant, e.g. in a wind. 

The minimum rotation rate as a fraction of breakup 



velocity A n 



[Sini/uc-b can also be checked using the 



same estimates for radius (Fig. llH . All objects are currently 
slow rotators (u eq sin i for FQ Aqr is probably overestimated 
for the reasons given above). 



Figure 11. Projected angular rotation as a fraction of breakup 
velocity A as a function of relative inverse radius f3 for extreme 
helium stars. 



In other words, the observed rotation rates for extreme 
helium stars are consistent with homologous contraction of 
the envelope at constant luminosity. They are not consistent 
with complete conservation of angular momentum during 
the white dwarf merger. With angular velocities less than 
one third the predicted rates, the likelihood of breakup as 
these stars contract towards the white dwarf phase is much 
reduced compared with Fig. |5| 
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10 CONCLUSION 

Assuming that angular momentum is strictly conserved 
within the system, we have calculated the rotation veloci- 
ties of stars produced by the merger of white dwarf binaries 
over a range of interesting initial mass ratios and possible 
outcomes. These include evolution through the giant phase 
and subsequently towards either the white dwarf or the he- 
lium main-sequence. We have demonstrated that while it 
may be possible to produce the initial giants, the star (or its 
core) will spin up as material is processed into the core or 
as the overall star contracts towards either of the compact 
configurations. This spinup would be high enough to cause 
breakup. 

More critically, the rotation velocities predicted under 
total conservation of angular momentum for post-giant ob- 
jects are greater than those observed in the putative CO+He 
merger products, the extreme helium stars, by a factor > 3. 

Therefore, for the merger model to work, the star must 
lose angular momentum at some point in its evolution. The 
most likely time for this to occur is following formation of 
the disk after breakup of the secondary. In order to assume 
total conservation of angular momentum, it was necessary 
to assume that this disk is collisionless, but that assump- 
tion would also prevent the accretion of material onto the 
primary. 

In fact the completely opposi t e ass umption must be 
considered. iLvnden-Bell fe Pringld (1197411 argue that, in a 
viscous disk, angular momentum will be transported out- 
wards while mass is transported inwards, on the viscous 
timescale T v j sc w R 2 /v, where R is the disk radius and v 
the viscosity. The latter is approximately ^c s H, where c s is 
the sound speed and the scale height H ~ Rc s /v c , v c being 
the Keplerian velocity for a circular orbit at radius R having 
a period Pr. Rearranging, we obtain r v i sc ~ 3R/v c (R/H) 2 . 
If we take the disk radius to be approximately four times 
the white dwarf radius R ~ 4J? w( j, and the disk scale height 
to be twice the white dwarf radius, we obtain v c « 6-R/Pr 
and Tvisc ~ 2Pr ~ 500s. Such a mechanism provides a very 
efficient sink for the angular momentum, but poses a prob- 
lem for the mass, which will be dumped onto the surface 
of the white dwarf at a rate far in excess of the Eddington 
rate (MEdd)- While accretion rates significantly higher than 
M^dd may be accommodated in a nonspherical geometry, 
e.g. by accretion at the equator balanced by radiation in 
the polar axis, the disparity here seems to be insurmount- 
able. Mass must be stored somewhere in the system where 
it is not supported by hydrostatic forces. 

A hybrid s olution might derive from a m echanism sug- 
gested fir st by [Lynden-Bell fe Pringld lll974h and explored 
again by |PophanifcNaravanrill99 ll) . This considers how 
material from an accretion disk would spin up the accre- 
tor (note, this is a true accretion disk, and not a disrupted 
secondary). When the accretor reaches critical velocity, it 
interacts with the disk, actually spinning up the disk. If this 
were to happen, then we could get momentum back into the 
disk, propagate it outwards and eject it from the outer edge 
of the disk. 

A more detailed study of the dynamics of the 
binary disruption, disk formation and subsequent ac- 
cretion is therefore warranted. Some progress has 
been made with smoothed particle hydrodynamics 



jBenz et alJ ll99fJ: |Se^retaui J _Clmbrierfc Mochkovitchl 
Il997t iGeurrero. Garcia- Be^^^senTT2o"o4l) ; calculations 
have dealt principally with the dynamics of the secondary 
disintegration and subsequent disk formation, and give 
useful insight into the disk heating and mass distribution, 
i.e. no explosion and negligible mass loss from the system. 
However, none goes far enough to establish the long-term 
outcome during and following disk accretion. 



ACKNOWLEDGEMENTS 

K. Gourgouliatos is grateful to the Armagh Observatory for 
a summer studentship during which this work was carried 
out. Research at the Armagh Observatory is funded by the 
Northern Ireland Department of Culture, Arts and Leisure. 
The authors are grateful to Profs D. Lynden-Bell and M. E. 
Bailey for stimulating discussions. 



REFERENCES 

Benz, W., Bowers, R. L., Cameron, A. G. W., Press, W. H. 

1990, ApJ 348, 647 
Chandrasekhar S. 1958, An Introduction to the Study of 

Stellar Structure. Dover Publications, New York, NY 
Chau W. Y. 1978, ApJ 219, 1038 
Eggleton P. P. 1983, ApJ 268, 368 

Geurrero, J., Garci'a-Berro, E., Isren, J. 2004, A&A 413, 
257 

Han, Z., Webbink, R. F., 1999, A&A 349, L17 
Iben I.,Jr. 1990, ApJ 353, 215 

Iben I., Jr., Tutukov A. V., Yungelson L. R., 1996, ApJ 456, 
750 

James R. A. 1964, ApJ 140, 552 

Jeffery C. S. 1993, A&A 279, 188 

Jeffery C. S., Heber U., 1993, A&A 270, 167 

Jeffery C. S., Hamill P. J., Harrison P. M, Jeffers, S. V. 

1998, A&A 340, 476 
Landau L., Lifshitz E., 1958, The Classical Theory of 

Fields. Pergamon, Oxford 
Lynden-Bell D., Pringle J. E. 1974, MNRAS 168, 603 
Marsh T. R. 1995, MNRAS 275, 1 
Mohan O, Bayaty A. R. 1980, Ap&SS 73, 227 
Napiwotzki, R., Karl, C. A., Nelemans, G., et al., 2005, In 

14th European Workshop on White Dwarfs, ASPC 334, 

375, eds. D.Koester, S.Moehler 
Nelemans O, Yungelson L. R., Portegies Zwart S. F., Ver- 

bunt F., 2001, A&A 365, 491 
Paczynski B. 1971, Act. Astr. 21, 1 
Paczynski B. 1990, ApJ 365, 9 

Pandey G., Lambert, D. L., Jeffery, C. S., Rao, N. K. 2006, 
ApJ 638, 454 

Panei J. A., Althaus L. O, Benvenuto O. G. 2000, A&A 
353, 970 

Popham R., Narayan R. 1991, ApJ 370, 604 
Pringle J. E. 1981, ARA&A 19, 137 
Saio H., Jeffery C. S. 1988, ApJ 328, 714 
Saio H., Jeffery C. S. 2000, MNRAS 313, 671 
Saio H., Jeffery C. S. 2002, MNRAS 333, 121 
Segretain, L., Chabrier, G., & Mochkovitch, R. 1997, ApJ 
481, 355 



On the angular momentum evolution of merged white dwarfs 11 



Tassoul J. L. 1978, Theory of Rotating Stars. Princeton 

University Press, Princeton, New Jersey 
Vennes S., Fontaine G., Brassard P. 1995, A&A 296, 117 
Webbink, R. F. 1984, ApJ 277, 355 



